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1 .  Introduction 

Using  a  local  balance  law  developed  by  Gurtin  [5]  for  the 
cohesive  zone,  we  give  a  system  of  equations  governing  the 
temperature  rise  around  a  fast  running  crack.  We  note  that  the 
effect  of  surface  temperature  cannot  be  neglected  for  a  fast 
running  crack,  and  we  discuss  tl. •:  equations  under  simplifying 
assumptions.  In  particular,  assuming  that  the  surface  temperature  is 
much  higher  than  the  surrounding  temperature,  we  arrive  at  a  simple 
solution  in  closed  form.  This  solution  agrees  with  experimental 
results  of  Puller,  Fox,  and  Field  [2]  for  polymethyl  methacrylate 
showing  that  the  temperature  rise  at  the  crack  tip  is  independent 
of  crack  speed. 
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2 .  Theory 


Consideration  will  be  restricted  to  a  semi- infinite  crack 
in  an  infinite  plate  stressed  symmetrically  with  respect  to 
the  plane  of  the  crack. 


We  assume  that  the  crack 
is  moving  with  constant 
velocity  v  in  the 
direction  of  negative  x. 
(See  the  Figure.) 


Neglecting  thermo-mechanical  coupling  we  have  the  following  ener 
balance  law  governing  the  temperature  rise  0  -  0(x,y,t)  away 


from  the  crack: 


i!  _  v t*2!  ■  a2e. 

ax2  ay2 


(2.1) 


where,  p  is  the  density,  k  conductivity  and  c  specific  heat. 

Let  z(t)  label  the  end  of  cohesive  zone  and  t  its  length,  so 
that  z(t)  +  l  gives  the  position  of  the  crack  tip. 

Since  the  work  rate  f(x,t)  due  to  plastic  work  in  the 
cohesive  zone  will  be  symmetric  about  the  plane  of  the  crack,  there 
will  be  no  heat  conduction  across  the  line  {y  *=  0,  x  <  z(t)); 


hence 


*  M  _  0 

*  ay  "  °' 


y  “  0,  x  <  z (t) . 


(2.2) 


We  allow  the  crack  surface  to  have  (surface)  temperature  rise 


(x,t) (>  C),  continuous  in  x(>  z(t)). 


According  to  a  recent  theory  of  Gurtin  [5] ,  the 


local  balance  law  in  the  cohesive  region  is  given  in  the  form 

ef  =  h  +  ,  (2.3) 

where  h  is  the  heat  flow  per  unit  surface  into  the  crack 
surface  from  the  body  (ef.,  e.g.,  (5.6)  in  [5]). 

Using  (2.3)  we  have  (by  routine  assumptions  and 
derivations) 


|  o(x,  t)  »  h  +  y  -  0,  x  >  z(t) , 

k  “  h»  y  =  0*,  x  >  z(t) , 


(2.4) 


where  p  is  a  constant.  We  assume  that  h  has  the  form 

h  =  a(0-«p) ,  (2.5) 

where  a  (*=  constant)  is  the  surface  conductivity.  We  note 
that,  in  (2.4)2  f(x,t)  =  0  at  x  >  z(t)  +  K,  and  we  neglect 
heat  loss  into  the  surrounding  air.1  In  Gurtin  [5],  the  energy 
of  the  newly  formed  free  surfaces  is  assumed  to  be  constant 
when  these  surfaces  are  exposed  to  a  constant  environmental 
temperature.  For  a  fast  running  crack,  however,  this  assumption 
is  redundant. 

Assuming  steady- state  conditions  and  introducing  the  moving 
coordinate  5  =  x  -  z(t),  the  equations  (2.1)- (2.5)  become 


XSee  Doll  [4]  . 
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pev  ||  «  k{^-§  +  ,  y  >  0,  -cd  <  S  <  oq 

dy 

k  ||  s  0,  y  =  0,  ?  <  0, 

|  v  If  -  a(0-o)  +  y  =  0,  ?  >  0, 


(2.6) 


*  f|  ■  a(0-cp) , 


y  -  o  ,  z  >  o. 


where  g(S)  s  *=  f(5  +  z(t),t)  is  the  work  rate  in  the  cohesive 


region.  By  (2.6)^  we  have 

b  =  30  +  RBI  (2.7) 

°  “  ay  +  2k  »  1  *  ; 

where  b  =  .  We  note  0  is  the  heat  required  to  raise 

2  o 

the  temperature  of  a  unit  surface  by  one  degree  (cal/fcm  /  C)  . 

We  expect  0  to  be  negligible  small.  The  dimensionless  constant 


b,  however,  is  approximately 

6  .  lx 10^X0  -  2024  Aluminium  Alloy 
1.7XlO5X0  -  Mild  Steel 
1.4XlO®X0  -  6A1-4V  Titanium  Alloy 

5.OX1O®X0  -  Polymethyl  Methacrylate  (k=5xl0-^/fcn/°C/s) 

wc  O 

for  v  =  500  m/sec.  Thus,  even  for  0  as  small  as  10  *-10  cal/ 

cmV°C,  we  may  not  neglect  the  term  b  |—  compared  to  in 

(2.7)  if  the  crack  velocity  is  as  large  as  500  n/s. 


We  will  discuss  the  temperature  rise  around  a  fast  running  crack 
using  the  above  equations,  but  under  certain  simplifying  assumptions . 
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3.  Analysis  and  Discussions 

Here  we  assume  that  6  «  co  on  the  crack  surface.  (Experi¬ 
ments  on  polymethyl  methacrylate  show  the  temperature  rise  on 
the  crack  faces  («l^m  in  depth)  to  be  about  500  K  throughout 
the  velocity  range  200~650  m/s  [2] ,  in  contrast,  the  maximum 
temperature  rise  at  a  distance  of  about  0.2~lmm  from  the  crack 
path,  is  only  0.1~1  K  [3]  .) 

Under  this  assumption  the  equations  (2.7)  become 


30 


^  y  i  » u  sj 

<>CVS?“  k(— 5 


2  2 
30  . 


ar  3y 


)  ,  y  >  0,  -0D  <  5  <  CD, 


y  «  o,  e  <  o 


I  v  If  =  -o»  +  ap,  y  =  0,  S  >  0, 


.  38 
k 


y  =  o*,  ?  >  o. 


We  assume  that  0  and  tend  to  zero  at  infinity.  We  then 


have  solutions  in  the  form 


(  ft 


(5-t) 


o 

.5 


2a 

ait).  ~  sv 
Bv  e 


-  P-t) 

SLltl  e  Pv 

Bv  e 


dt  t  £  5, 

dt  0  £  ?  < 

5  <  0, 


where 


0(5, y)  =  1°°  K(5-t,y) dt. 


\  1  ax  ,  /  2,  2.  pev 

K(x,y)  =  -  e  Kc(a  >'x  +y  ) ,  a  = 


i 


4 

I 
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with  K  the  modified  Bessel  function  of  the  second  kind, 
o 

Using  the  Dugdale  model,  the  work  rate  due  to  plastic  work 
is  given  in  the  form* 


where  o_ 
o 

Thus  the 


and  thus. 


-  °o 


is  the  yield  stress  and  6  the  separation  distance, 
temperature  rise  «o{?)  for  %,  ^  l  becomes 


«p(5)  **  (e 


2al 

vp 


I 


.  _  d5(t)  2a  ^  2a,  „  .. 


o  .dt 

p 


,vp 


dt)  e 


vp 


2fl*  ,  _  d6  (jop  208... 


**  (e 


vp 


dx 


P 


2a 

vp  ”,  %  "  vp 
e  K  dx) e  v 


(5~t) 


(3.1) 


if  2a4,/vp  is  negligibly  small,  (3.1)  is  approximately, 

oj  (*>  -  ^ 


*(?> 


Pv 


(3.2) 


where  r?  =  ?-<(.. 


This  result  shows,  interestingly,  that  the  temperature  rise 
at  the  crack  tip  (r>  =  0)  _is  proportional  to  rrQ6  (t)  with 
constant  of  proportionality  is  l/p . 


*The  work  rate  g(£)  for  the  Dugdale  model 
is  computed  by  Levy  and  Rice  [1)  as 


g(?)  = 


4(l-v2)o2v 


Ln(l  ■+ 

1  -  J?7l 


in  small  plane  strain 


0  <  t  <  l. 


2 

Note  that  we  have  assumed  plastic  work  in  the  cohesive 
zone  is  completely  converted  into  heat.  (This  is  gener¬ 
ally  true  of  the  energy  expended  in  plastically  deforming 
a  metal.)  For  partial  conversion  g  should  be  scaled 
down  appr  opr i at e ly . 
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Note  that  tp  depends  on  v  only  through  a  dependence 
on  u/v.  For  the  point  x  =  y  =  0,  say,  occupied  by  the  tip  at 
t  —  0,  t)  will  be  the  distance  frcm  the  tip  at  time  t  =  t)/v. 

Thus,  by  (3.2),  the  temperature  rise  at  this  point  should  depend 
only  on  time;  it  should  be  independent1  of  v.  This  result  is  in 
agreement  with  experimental  result  of  Fuller,  Fox,  and  Field  [2], 
who  found  that  the  temperature  rise  at  a  fixed  point  on  the  axis 
of  the  crack  was  approximately  independent  of  crack  velocity. 

In  fact,  they  note  that  "the  results  combined  to  give  a  temperature 
rise  of  approximately  500  K  throughout  the  velocity  range  studied 
(200-650  xo/s)  ." 


For  the  fixed  point  x  =  y  =  0,  (3.2)  gives 

Ob  (t)  -  -^t 

<P  «  _o_ -  e  8 


and  hence  «p  has  the  form 


*  -Bt 
tp  w  Ae  , 


where  A  =  oq6 (l)/P  and  B  =  -2a/8. 

The  data  of  [2] ,  when  averaged,  give 


tp  =  457  K 

at 

t  =  10us 

tp  =  361  K 

at 

t  =  20us 

tp  =  304  K 

at 

t  =  35|is. 

Using  the  values  at  t  =  10us 

and 

t  =  35us 
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A  =  5.38x102k, 

B  =  1.63Xl04sec 


we  find  that 


Here  we  follow  Levy  and  Rice  [1)  and  assume  that 
aQ6 (t)  is  independent  of  v.  Levy  and  Rice  [1] ,  however, 
found  that  this  temperature  rise  is  proportional  to  /v. 
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This  gives  <p  =  388  K  at  t  =  20us  as  compared  to 
the  value  to  =  361  K  of  [2]  .  Further,  we  find  that  cp  *=  538  K 
at  t  =  C;  that  is,  the  temperature  rise  at  the  crack  tip  is 
538  K.  For  the  values  l  -  1  mm  (which  we  feel  is  an  upper 
bound  for  t)  and  v  =  200  ir/s,  we  find  that  the  dimensionless 
constant 

2at/vP  =  tB/v, 

which  we  neglected  in  defining  (3.2),  has  the  approximate  value 
of  0.08. 
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